A single-spin precessing gravitational wave in closed form 
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Current searches for gravitational waves from compact binaries use templates where the binary 
does not precess. By contrast, generic black hole-neutron star binaries will significantly precess, 
inducing modulations that nonprecessing searches are particularly ill-suited to recover. In this 
paper we provide a closed-form representation of the single-spin precessing waveform in the frequency 
domain by reorganizing the signal as a sum over harmonics, each of which resembles a nonprecessing 
waveform. This form enables simple analytic calculations (e.g., a Fisher matrix) and computationally 
fast frequency-domain templates. We have verified that our SpinTaylorF2SingleSpin model is 
faithful to the time-domain SpinTaylorT2 model: for the majority of generic BH-NS binaries, the 
corresponding signals generated from each model agree to better than 1%. Our waveform model is 
now available as part of the publicly-available lalsimulation waveform-generation code. 

PACS numbers: 



For the astrophysically most plausible strong gravi- 
tational wave sources - coalescing compact binaries of 
black holes and neutron stars with total mass M — 
mi -\- TO2 < I6M0 - ground based gravitational wave 
detector networks (notably LIGO [1] and Virgo |2l E]) 
are principally sensitive to each binary's nearly-adiabatic 
and quasi-circular inspiral 4-12]. For nonprecessing bi- 
naries, the inspiral signal model has been explored in 
detail, both analytically and with numerical simulation. 
Fast, stationary-phase approximations to the signal exist 
which faithfully reproduce the signal [T3]. By contrast, 
plausible astrophysical processes will produce merging 
black hole binaries with arbitrary spin orientations, fully 
populating the 15-dimensional model space for quasi- 
circular inspiral. Fully generic binaries undergo ex- 
tremely complicated evolution, consisting of a gradually 
decaying and rapidly precessing orbit |14H16j . For many 
astrophysically interesting sources, including many bina- 
ries containing at least one black hole, significant spin 
precession is expected while the signal passes through 
LIGO and Virgo's sensitive frequency band. Precessing 
sources have been substantially less-thoroughly modeled. 
At present, the current state-of-the-art SpinTaylorT4 
O [6] and SpinTaylorT2 models [13] solve a coupled 
time-domain ordinary differential equation to evolve the 
orbits and spins. Until this work, no fast, stationary- 
phase approximation existed which faithfully reproduces 
a precessing gravitational waveform. The need to nu- 
merically evolve several differential equations imposes a 
significant computational cost. The added physical pa- 
rameters needed to capture precessing spins multiply the 
computational burden of using precessing signals for data 
analysis. On the one hand, search strategies attempting 
to recover precessing signals could not employ the full sig- 
nal, instead adopting a detection template family, whose 
waveforms usually at best imperfectly approximate the 
target waveforms O [7J |H1 UHl E] but whose template 
banks are smaller and better understood. On the other 
hand, the large signal space required an equally large 
number of slow-to-evaluate templates to cover it, placing 



a similarly-large computational burden on gravitational 
wave parameter estimation using physical signals. Fi- 
nally, the Fisher matrix and other quantities arising nat- 
urally in data analysis take a simpler form in our model, 
where derivatives can be efficiently evaluated [18]. By 
contrast, time-domain models are essentially opaque to 
analysis, since interesting calculations can only be done 
numerically. 

In this paper, we show that the gravitational wave sig- 
nal for a single-spin binary can be well approximated with 
a simple analytic form in the frequency domain. Our de- 
composition expresses the signal as a sum of five terms, 
each a weak modulation (sideband) of a leading-order 
( "carrier" ) function: 
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Critically, our decomposition provides closed-form func- 
tions of frequency for all of the terms entering this ex- 
pression. £",„(/) can be thought of as a complex antenna 
pattern which contains all of the orientation dependence 
of the waveform, and evolves very slowly over radiation 
reaction rather than orbital timescales. 

We have tested our frequency-domain 
SpinTaylorF2SingleSpin signal model to the time- 
domain SpinTaylorT2 approximant J13j by calculating 
the faithfulness. The noise-weighted inner product of 
two waveforms hi and /i2 is [19J 
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and the overlap is {hi\h2) / y/ {hi\hi){h2\h2) . We gener- 
ate both signals with the exact same mass, spin, and 
orientation parameters, then maximize the overlap over 
only the time and phase, yielding the so-called faithful- 
ness. The same polarization of the two waveforms is 
compared. This is a more stringent requirement than 
requiring that our model be effectual (high overlap when 
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FIG. 1: Faithfulness of our model; For 10 + IAMq 
BH-NS binaries with random orientations and random BH 
spins with magnitudes in [0.5,1], we compute the faith- 
fulness (overlap maximized over time and phase) of our 
SpinTaylorF2SingleSpin model with SpinTaylorT2. The 
agreement is better that 2% except for a small number of 
cases with very large spin-orbit misalignment. 



maximizing over parameters like masses, as in a search). 
We did 5000 simulations of 10 + IAMq BH-NS binaries 
with random source orientations, and BH spin with ran- 
dom direction and magnitudes x S [0.5,1]. For 5'„(/), 
we use the Advanced LIGO high-power zero-detune noise 
spectrum [Ml 1^ with /; = 15Hz and the ISCO frequency 
as fh- The two models agree to better than 2% in almost 
all cases, with mismatches of up to 10% occuring only 
for strong spin-orbit misalignment {L ■ S — k < —0.9). 
This demonstrates that the approximations that went 
into building our model are all justified. 

In our notation, the components of the binary have 
masses mi and m2, with mi the larger mass; for conve- 
nience, we also use M = mi + m2, rj = mim2/M^, and 
the mass ratio q = 1112/ mi. Only the larger mass has 
spin, with a dimcnsionless spin parameter xi s-'^d spin 
I Si I = mfxi- The orbital angular velocity is denoted 
Lj and is related to the gravitational wave frequency by 
uj = irf and to the orbital velocity by u = {Mu))^/'^. 

Because the waveform only has two independent po- 
larizations along any line of sight, we can without loss of 
generality treat it as overhead of the detector. Any other 
sky position just acts as an amplitude scaling and a rota- 
tion of the polarization. To define the orientation of the 
binary with respect to the detector, we use the total an- 
gular momentum J = L -|- S which is nearly fixed during 
the inspiral. The inclination of J to the line toward the 
detector is 0, and the polarization angle -0 is defined by 
the angle between J and the x-arm of the detector in the 
plane of the sky; see [22] , [13] , and references therein. 

The precessing waveform can be decomposed into a 
non-precessing waveform acted on by a time-dependent 
rotation [M1 - E7] . The rotation accounts for the precession 
of the orbital angular momentum of the binary. The 



waveform in the source frame is 
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where $ is the orbital phase and where here and hence- 
forth ~ denotes the frame co-rotating with the precession 
of L. In this frame, the source amplitudes nearly obey 
hi.-m — i^^yi^em-' ^'^ ^^^ extent spin-dependent higher 
harmonics can be neglected; expressions for /i'™ exist in 
the literature [1]. 

The source frame and inertial frame are related by a 
time-dependent rotation, specified by three Euler angles 
(a, /3, 7).^ The source frame is aligned with the instanta- 
neous angular momentum L; the inertial frame is aligned 
with the total angular momentum J of the binary. In ad- 
dition to the two Euler angles (a, /3) needed to map J — > 
L, we apply a third Euler angle 7 = — / cos /3 {da/ dt)dt to 
minimize superfluous coordinate changes associated with 
the orbital plane; see Boyle et al. [23]. By applying 
this time-dependent rotation, the waveform in an iner- 
tial frame can be expressed as a weighted sum of the co- 
rotating- frame amplitudes /i^, this rotation allows can 
be decomposed as the waveform in a frame co-rotating 
with the precession and a time-dependent transformation 
from the co-rotating to the detector coordinates. 
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where D^',^ is the usual Wigner rotation matrix rep- 
resentation of SU(2) ^24'-'26]. Because ft.+ — ihx has 
spin weight —2, this expression must be proportional to 
exp(—2iTp), where the polarization angle -0 is the angle 
of J projected into the plane of the sky. 

For the special case of leading-order quadrupole emis- 
sion (/i^™ = unless I — \m\ = 2), the real part of the 
above expression can be rewritten as [22j 
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where z{q) is a complex variable that captures the mod- 
ulation due to precession. These effects are often split 
into amplitude and phase modulations [14!, or the wave- 
form is decomposed in spherical harmonics on 9 and ^ 



^ The three Euler angles (a,/3,7) can be visualized as follows. 
Start with the source frame aligned with the detector frame, 
so that L is along the z axis of the inertial frame. Rotate by 7 
around z, which physically contributes only to the orbital phase. 
Then rotate by fi around y, which inclines L relative to z. The 
final rotation is by a around z, which accounts for the preces- 
sion of L around J. Our angle sign convention follows |24| and 
disagrees with I28| . who use an adjoint representation. 



[5]. Instead, it is possible to write 
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where the latter expression implicitly defines the coeffi- 
cients Zm- The Zm are normalized so Zm = (5m. 2 when 
Q = ip = 0, In this sum, each term is proportional to 
giiTiQ g^^^ jg therefore modulated at a harmonic of the pre- 
cession frequency; multiplying the leading-order exp2i$, 
each term therefore produces a sideband, offset from the 
carrier frequency. In general the i?,„ are given by 
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While this expression applies in general, we substantially 
increase its value by adopting coordinates with a hierar- 
chy of timescales, so $ changes on an orbital timescale; 
a on a radiation reaction timescale; and /3 and v on the 
inspiral timescale. This separation of timescales occurs 
naturally in single-spin BH-NS binaries. To an excel- 
lent approximation, a single-spin binary undergoes sim- 
ple precession |14j . where the total angular momentum 
direction is fixed; on short (precession) timescales the or- 
bital angular momentum precesses around J at a uniform 
rate; and on longer (inspiral) timescales the precession 
cone opening angle gradually increases. Adopting coor- 
dinates aligned with J, the polar angle /3 is identified 
as the precession cone opening angle [55] and therefore 
changes slowly. 

In the co-rotating frame, orbital evolution can be cal- 
culated using the instantaneous binary's binding energy 
E{v) and the gravitational wave flux J-{v) |BJ: 
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Both E{v) and T{v) are known as post-Newtonian ex- 
pansions in the velocity v. These are currently known to 
order v^ past the leading order for non-spinning terms 
and v^ in terms involving spin (note that terms at order 
v^ and beyond may also have terms containing powers 
of logii). Since dt/dv is already approximate, we are not 
unjustified in expanding dt/dv as a power series in w, to 
the same order as we know the flux and energy, to flnd 
t{v) as a closed-form series in v. We use the additional 
relation for the orbital phase 
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to obtain $(w) in closed form. The gravitational- wave 
frequency is related to the velocity by v = (ttM/)^/^ by 
Kepler's Law and the fact that the gravitational-wave 
frequency is twice the orbital frequency. 



For the special case of BH-NS binaries, we can analyt- 
ically solve the spin-precession equations that determine 
L and hence a, /?, and Q as functions of the velocity. 
These three angles can be substituted into the rotation 
operator to transform the waveform from the precessing 
source frame into the inertial detector frame. 

The leading order ("Newtonian") expression for the 
orbital angular momentum is 
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and for the total angular momentum is 
J = L + Si . 



(10) 



(11) 



The magnitude of the spin is conserved. In terms of the 
dimensionless spin, it is Si = mlx^- The angle between 
L and Si does not evolve in the single-spin case, giving 
the conserved quantity k = Ljv 'Si. We can additionally 
define two ratios 
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The global behavior of Tj is not well-fit by a single low- 
order polynomial, as it approaches 1 at w — >■ and is 
proportional to v in the limit of large v. As a result, 
expressions involving T j are generally not well-fit by a 
standard PN expansion. 

In terms of these quantities, the opening angle /3 of the 
cone swept out by L (the "precession cone") is 
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The time evolution of La? is given by 
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which causes Lat to rotate around J with an angular rate 
of 



The definitions of a and Q are then 

C, = a cos (3 — Vtp cos /? . 



(16) 
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Substituting the definitions used above leads to a simple 
expression 



4 



We use the TaylorT2 [5] expression for dt/dv as a power 
series in v to express a and (^ as integrals in v rather 
than in t. Despite the non-polynomial behavior in Tj, 
both integrals can be evaluated term-by-term to produce 
closed-form expressions for a{v) and C{v) at least to 3 PN 
order, where terms of the form u" log v appear. These 
closed-form expressions are provided as supplementary 
online material, and in the lalsimulation code. Since 
a, C and the orbital phase [Eq. ^] influence each terms' 
phase in a similar way but a, C depend on integrals over 
dt/dv multiplied by at least v^ , at least two fewer terms in 
a post-Newtonian series expansion for dt/dv are needed 
to reproduce the precessional dynamics of a, (^ at the ac- 
curacy needed. 

For nonprecessing signals, a commonly-used ap- 
proximation to a Fourier transform is provided 
by the stationary-phase approximation, which for 
a real-valued signal ReAexp(i$(i)) has the form 
A{v{f)) /2yJ id^f^ / dt"^ exp #(/)> where where V(/) = 
27r/i(v) - 2$(z;); see, e.g., Eq. (2.18) in [19]. Using 
this approximation, the Fourier transform of our mod- 
ulated waveform [Eq. ([5])] can be efficiently computed 
term-by-term, using the phase ^m = 2(<f> — 7) -I- ma. 
Following custom in gravitational wave data analysis, 
we approximate the stationary phase amplitude factor 
ex l/y/(PWJdt/^ ~ (d2$/di2)-i/2 Yyy its leading order 
term, proportional to /^^/^. We therefore find 
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where the expressions for a{v{f)) and "f{v{f)) follow by 
explicit substitution. 

Our fast, faithful frequency-domain waveform pro- 
vides an efficient, easily-understood, and analytically- 
tractable way to represent the gravitational wave signal 
from generic precessing BH-NS binaries. Our method 



generalizes naturally to include higher harmonics. As 
noted above, this method allows for fast, accurate an- 
alytic Fisher matrices, to estimate the performance of 
parameter estimation. More broadly, overlaps between 
generic single-spin sources can both be calculated and 
understood analytically. As a concrete example, Brown 
et al. P5j evaluated the overlap between a generic pre- 
cessing signal and a nonprecessing search template. In 
our representation, each sideband has a unique time- 
frequency trajectory, offset from the "chirp" of the or- 
bital frequency versus time. Hence, the best fit between 
a nonprecessing model will occur when the nonprecess- 
ing model's time-frequency path lies on one of the side- 
bands; a nonprecessing search misses all power, except 
that associated with the optimal sideband; and the best- 
fitting match will be proportional to \zm\, reproducing 
their general result. Finally, our model suggests several 
strategies for a viable search for precessing single-spin 
signals. First and foremost, our model allows us to iden- 
tify precisely which masses, spins, and viewing orienta- 
tions would benefit from a multi-modal search. Second, 
the functional form of z^ ensures that rarely more than 
two m will contribute significantly to the signal. Each 
of the harmonics is essentially a nonprecessing waveform 
and will likely match well an aligned-spin waveform with 
different parameters. An analytic aligned spin template 
bank was presented in |29| . Extending the boundaries 
of the bank would allow it to pick up each of the har- 
monics as an isolated spot of power. A promising search 
strategy would be to recombine the power from two (or 
more) spots of power and so recover most of the power in 
the precessing signal. Third, our expression simplifies the 
form of and overlaps between the I, m modes of precessing 
signals. Combined with a massive reduction in compu- 
tational cost, our model may allow effective searches for 
generic precessing systems, using physical templates and 
maximizing over source orientation [5]. 
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